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ABSTRACT 

AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU-A and HSB, to form a next generation polar 
orbiting infrared and microwave atmospheric sounding system. AIRS is a grating spectrometer with a number of linear 
arrays of detectors with each detector sensitive to outgoing radiation in a characteristic frequency v t with a spectral band 
pass Av* of roughly Vf/1200. AIRS contains 2378 spectral channels covering portions of the spectral region 650 cm' 1 
(15.38 (im) - 2665 cm' 1 (3.752 p,m). These spectral regions contain significant absorption features from two C0 2 
absorption bands, the 15 |im (longwave) C0 2 band, and the 4.3 pm (shortwave) C0 2 absorption band. There are also two 
atmospheric window regions, the 12 pm - 8 jim (longwave) window, and the 4.17 jim - 3.75 pm (shortwave) window. 
Historically, determination of surface and atmospheric temperatures from satellite observations was performed using 
primarily observations in the longwave window and C0 2 absorption regions. One reason for this was concerns about the 
effects, during the day, of reflected sunlight and non-Local Thermodynamic Equilibrium (non-LTE) on the observed 
radiances in the shortwave portion of the spectrum. According to cloud clearing theory, more accurate soundings of both 
surface skin and atmospheric temperatures can be obtained under partial cloud cover conditions if one uses the longwave 
channels to determine cloud cleared radiances for all channels, and uses only from shortwave channels in the 
determination of surface and atmospheric temperatures. This procedure is now being used by the AIRS Science Team in 
preparation for the AIRS Version 6 Retrieval Algorithm. This paper describes how the effects on the radiances of solar 
radiation reflected by clouds and the Earth’s surface, and also of non-LTE, are accounted for in the analysis of the data. 
Results are presented for both daytime and nighttime conditions showing improved surface and atmospheric soundings 
under partial cloud cover resulted from not using R t in the retrieval process for any longwave channels sensitive to cloud 
effects. This improvement is made possible because AIRS NEDT in the shortwave portion of the spectrum is extremely 
low. 

Keywords: AIRS/AMSU, high spectral resolution IR sounders, retrieval methodology, IR sounding in cloudy 
conditions, cloud cleared radiances. Quality Control. 

1. INTRODUCTION 

AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU-A and HSB, to form a next generation polar 
orbiting infrared and microwave atmospheric sounding system 1 . AIRS is a grating spectrometer with a number of linear 
arrays of detectors with each detector sensitive to outgoing radiation in a characteristic frequency with a spectral band 
pass AVf of roughly Vj/1200. AIRS contains 2378 spectral channels covering portions of the spectral region 650 cm' 1 
(15.38 pm) - 2665 cm' 1 (3.752 pm). The spectral sampling interval (except for the existence of a few gaps) is vj 2400, 
giving two AIRS channels per spectral half width. AIRS is accompanied by the microwave instrument AMSU-A. There 
is a 3x3 array of AIRS footprints within a given AMSU-A footprint, with spatial resolutions of 13 km and 45 km at nadir 
viewing respectively. 

Figure 1 shows NEDT in blue for AIRS channels evaluated at the channel brightness temperatures 0* for a clear sky 
tropical scene. Figure 1 also includes the noise in red for the IASI instrument evaluated at IASI 0j for the same scene. 
IASI was launched on October 19, 2006 on MetOp-A, also accompanied by AMSU-A 2 . IASI has 4 footprints within an 
AMSU-A footprint. IASI is an interferometer with a maximum Optical Path Displacement (OPD) of 2.0 cm. The spectral 
sampling of IASI is 0.25 cm" 1 throughout the spectrum. This is roughly the same as AIRS at 650 cm' 1 , and about 4 times 
more often than AIRS at 2400 cm' 1 . Channel noise for AIRS and IASI is comparable throughout most of the spectrum 
but IASI channel noise is very large relative to AIRS in the shortwave portion of the spectrum. In principle, better noise 


performance could have been obtained by IASI at these short wavelengths. Channels in this portion of the spectrum were 
not properly utilized by other researchers in the retrieval process because of concerns about effects on the observed 
radiances, during the day, of solar radiation reflected by both the surface and clouds, as well as effects of non-Local 
Thermodynamic Equilibrium (non-LTE). Therefore, noise performance in this portion of the spectrum was not a priority 
in the design of IASI. The methodology we use in this paper to determine improved surface and atmospheric temperature 
using only shortwave AIRS observations is not practical using IASI because the IASI noise performance is not adequate 
in the shortwave portion of the spectrum. 
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The sounding goals of AIRS are to produce 1 km tropospheric layer mean temperatures with an rms error of IK, and 1 
km layer precipitable water with an rms error of 20 percent, in cases with up to 80 percent effective cloud cover. The 
primary products of AIRS/AMSU-A are twice daily global fields of atmospheric temperature-humidity profiles, ozone 
profiles, sea/land surface skin temperature, and cloud related parameters including OLR. Also included are the 
reconstructed clear column radiances used to generate these products. is a derived quantity representative of the 
radiance AIRS would have seen in channel i if there were no clouds in the field of view. All products also have error 
estimates. The products are designed for data assimilation purposes so as to improve numerical weather prediction, as 
well as for the study of climate and meteorological processes. With regard to data assimilation, one can use either the 
products themselves or the clear column radiances from which the products were derived. From both the data 
assimilation and climate perspectives, better spatial coverage of the retrievals is desirable provided the accuracy of 
retrieved products is sufficient. 

The theoretical approach used to analyze AIRS/AMSU/HSB data in the presence of clouds in the AIRS Science Team 
Version 3 at-launch algorithm, and that used in the Version 4 post-launch algorithm, have been published previously 3,4 . 
Significant theoretical and practical improvements have been made in the analysis of AIRS/AMSU data since the 
Version 4 algorithm. Most of these have already been incorporated in the AIRS Science Team Version 5 algorithm 5 , now 
being used operationally at the Goddard DISC. The AIRS Version 5 retrieval algorithm contained three significant 
improvements over Version 4. Improved physics in Version 5 allowed for use of AIRS clear column radiances R t in the 
entire 4.3 pm C0 2 absorption band in the retrieval of temperature profiles T(p) during both day and night. Tropospheric 
sounding 15 pm C0 2 observations were used primarily in the generation of clear column radiances for all channels. 
This new approach allowed for the generation of accurate Quality Controlled values of fii and T(p) under more stressing 
cloud conditions. Secondly, Version 5 contained a new methodology to provide accurate case-by-case error estimates for 
retrieved geophysical parameters and for channel-by-channel clear column radiances. Thresholds of these error estimates 
were used in the new approach for Quality Control. Finally, Version 5 contained for the first time an approach to provide 
AIRS soundings in partially cloudy conditions that does not require use of any microwave data. This new AIRS Only 
sounding methodology, referred to as AIRS Version 5 AO, was developed as a backup to AIRS Version 5 should the 



AMSU-A instrument fail. Susskind et al. 5 show that Version 5 AIRS Only soundings are only slightly degraded from the 
AIRS/AMSU soundings, even at large fractional cloud cover. 

Significant further progress has been made by the AIRS Science Team since the delivery of the AIRS Version 5 retrieval 
algorithm in preparation of the AIRS Science Team Version 6 AIRS/AMSU retrieval algorithm expected to be made 
operational at the Goddard DISC in early 2010. The current version of the retrieval algorithm, as of July 2009, is referred 
to in this paper as Version 5.24. Version 5.24 is very similar to Version 5 with one major exception. As in Version 5, the 
generation of cloud cleared radiances R t for all channels is performed using observed radiances R t for longwave 15 Jim 
and 11 pm channels. Tropospheric temperature profiles were retrieved in Version 5 using only the AIRS shortwave 4.2 
|im C0 2 channels, but surface skin temperature was retrieved simultaneously with surface spectral emissivity and bi- 
directional reflectance using observations both in the longwave 8-12 pm window region and in the shortwave 4.0 Jim - 
3.76 |im window region. In Version 5.24, only window observations in the shortwave window region, 4.0 pm - 3.76 pm 
are used to determine surface skin temperatures as well as shortwave surface spectral emissivities and surface bi- 
directional reflectance. The current use of only shortwave AIRS channels in the retrieval of both atmospheric and surface 
parameters has resulted in significant improvement in the ability to obtain both accurate temperature profiles and surface 
skin temperatures under more stressing partial cloud cover conditions. This paper discusses some of the theoretical 
improvements made both in Version 5 as compared to Version 4, and Version 5.24 compared to Version 5. 


2. OVERVIEW OF THE AIRS RETRIEVAL METHODOLOGY 

The AIRS Science Team Version 5, Version 5 AO, and Version 5.24 retrieval algorithms are very similar to each other 
and to Versions of the AIRS Science Team retrieval algorithms described previously 3,4 . The basic approach and 
methodology used to determine geophysical parameters from observed AIRS/AMSU radiances is identical to that 
described in great detail in 3 . The basic approached is summarized below. 

Retrievals of all geophysical parameters are physically based and represent states Xj determined for case c that best 
match a set of clear column radiance R ic for the subset of AIRS channels i used in the retrieval process. The clear 
column radiances R i c are a derived quantity that represent the radiances that AIRS “would have seen” if the 3x3 AIRS 
field of regard (FOR) c were completely clear. Retrievals of geophysical parameters are performed sequentially, that is, 
only a subset of the geophysical parameters within the state Xj is modified from that of the incoming state A^in a given 
step. A GCM forecast is not used in any way in the retrieval procedure, except for the forecasted surface pressure p surf , 
which is used as the lower pressure boundary when computing expected radiances R i CL ^for a given geophysical state. 

In Version 5, the steps are as follows: After a start up procedure giving the state used to determine the initial clear 
column radiances AIRS/AMSU data are used to retrieve: a) surface skin temperature, surface spectral emissivity and 
surface bi-directional reflectance of solar radiation; b) atmospheric temperature profile; c) atmospheric moisture profile; 
d) atmospheric ozone profile; e) atmospheric CO profile; f) atmospheric CH 4 profile; and g) cloud properties. These steps 
are done sequentially, solving only for the variables to be determined in each retrieval step and using previously 
determined variables as fixed but with an appropriate uncertainty attached to them which is accounted for in the channel 
noise covariance matrix N. The objective in each step (a-f) is to find solutions which best match fii for the subset of 
channels selected for use in that step, bearing in mind the channel noise estimates. Steps a-f are ordered so as to allow for 
selection of channels in each step which are primarily sensitive to variables to be determined in that step or determined in 
a previous step, and relatively insensitive to other parameters. Separation of the problem in this manner also allows for 
the problem in each step to be made as linear as possible. Step g is performed after the surface and atmospheric 
conditions have been determined using the observed radiances R t . 

In Version 5.24, there is a slight modification to the sequence of steps used in Version 5, as there is a new step 
performed in the retrieval sequence. In Version 5, step a) used channels in both the longwave and shortwave window 
regions and simultaneously solved for surface skin temperature T s , shortwave surface spectral emissivity £ SV y(v) an< ^ 
surface spectral bi-directional reflectance p^ (v), and longwave surface spectral emissivity £{ W (v). In Version 5.24, 
only shortwave window channels are used in this retrieval step to simultaneously determine T s , £ sw (v), and p sw (v). The 


longwave surface spectral emissivity £^ w (v) is now solved for in a subsequent step using only channels in the longwave 
window spectral region. This new step is performed after the humidity profile retrieval step because longwave window 
radiances can be very sensitive to the amount of atmospheric water vapor. 


3. BASIC THEORETICAL CONSIDERATIONS FOR CLOUD CLEARING 


Using the assumption that adjacent fields of view have otherwise identical geophysical conditions except for cloud cover, 
Chahine 6 has shown that in the case of K-l cloud formations, observations in K fields of view (FOV’s) are needed to 
obtain channel i cloud cleared radiances fii which can be determined using the equation 

= Ri,l "J" 2/c= l Vk (^i,l — Ri,K+l-k) (1) 

In which R ik is the channel i radiance observation in field of view k. 

In analyzing AIRS/AMSU-A data, a single sounding is produced using all 9 AIRS FOV’s falling within a single AMSU- 
A footprint 3,4 . The AIRS Science Team has found it is advantageous to extrapolate the radiances in the 9 fields of view 
according to a similar equation of the form 


— Ri,AVG + £fc= 1 Vk ( Ri,AVG ~ R i,k ) 


( 2 ) 


where Ri tA vc i s the average channel i radiance of all 9 fields of view. The superscript n has been added to Equation 2 
because the generation of R t is iterative 5 . Optimal values of rj k will give true values of R t up to instrumental noise 
effects. While there are 9 values of rj k shown in Equation 9, only 8 of them are linearly independent. 

Most of the 7=59 cloud clearing channels used in Equation 2 are channels in the 15 jum C0 2 absorption band with 
temperature weighting functions covering the atmosphere from the lower stratosphere to the surface. Observations in 
some longwave window channels are used as well as the determination of rj k . 


3.1 Determination of rj 

If, for each channel i, one substitutes an estimate of for Rf in Equation 2, this gives 7 equations for K (=9) 

unknowns. The unconstrained weighted least square solution to this multilinear problem is given by 

Vkx 1 = (3) 

where A R is an IxK matrix with A R i k — R AVG — R i k , A R glr is an Ixl matrix given by A R* clr = R?,clr ~ R uavg> anc * A 
is an Ixl channel noise covariance matrix. R^clr ls generated by computing expected radiances for cloud clearing 
channel j based on the current estimate of the geophysical state X n assuming cloud free conditions. As in Susskind et 
al. 3,4 , the solution for 77 is stabilized by solving for coefficients of up to the first four principal components of 
[AR'N~ 1 AR]. The parameters rj k determined in Equation 3, using I channels, are assumed to characterize the cloud 
formations and thus should be valid for use in all channels. 

The key to the accurate determination of r| is obtaining the best estimates of AR” CLRi along with an accurate treatment 
of the noise covariance matrix N. The values of A R™ clr , which are used to determine r\ (and R t ) are iterative and are 
computed based on the current best estimate of all relevant surface and atmospheric properties. It is best for the estimated 
geophysical parameters to be unbiased over large regions of the atmosphere. In all previous versions of the cloud 
clearing algorithm, including Version 4, the geophysical state X° used to estimate Ri CLR was derived from an AMSU-A 
retrieval state, thus insuring an unbiased temperature and moisture profile over coarse layers in the atmosphere. 
Subsequent research has shown that a reasonably good regression relationship can be obtained between geophysical 
parameters and observed (un-cloud cleared) AIRS radiances, and this regression based state can be used to generate an 


alternative initial state X° used for initial cloud clearing. Generation of this state X° does not require use of any AMSU 
observations. This is the approach used to obtain X° in the AIRS Only cloud clearing system, Version 5 AO. Version 5 
uses analogous methodology to give the initial cloud clearing state, but the cloudy regression makes use of both AIRS 
and AMSU observations. The state derived from this cloudy regression, X reg is followed in Version 5 by an AMSU only 
temperature profile retrieval step to produce the state X° which is used for initial cloud clearing. Such a step is not 
performed in Version 5 AO. 

The cloudy regression can produce biased initial states, especially if it is not followed by an AMSU-A retrieval step. If 
the state 7°(p) used to derive R^clr were biased (say too warm), incorrect values of 7] k would be determined which 
would result in Rf being too large, which in turn would result in the retrieved 7’ 1 (p) being too warm. Chahine 7 has 
shown that it is optimal to use only longwave (15 Jim and 12 pm) channels for cloud clearing, and shortwave 4.2 pm 
channels for the determination of temperature profiles. This is done so as to minimize the bias in retrieved temperature 
profiles X n+1 resulting from biases in the temperature profile X n used to determine rj k . Up to Version 4, most 4.2 pm 
channels could not be utilized during the day because these channels are affected by non-Local Thermodynamic 
Equilibrium (non-LTE) which was not accounted for in the Version 4 Radiative Transfer Algorithm (RTA) 4 . Therefore, 
the optimal combination of channels for cloud clearing and sounding purposes was not used in Version 4, in which most 
temperature sounding channels were in the 15 pm CO 2 band. 

A significant improvement over the AIRS Science Team Version 4 retrieval algorithm is the use of a new RTA in 
Version 5 which contains improved characterization of atmospheric absorption characteristics, and more significantly, 
accounts for effects of non-LTE 8 . This for the first time enabled the use of almost all shortwave C0 2 channels in the 
temperature profile retrieval step of the physical retrieval algorithm, both day and night. This new RTA also contains a 
better parameterization of the absorption characteristics of the AIRS channels as a function of atmospheric temperature 
and constituent profile as a function of satellite zenith angle. Stratospheric 15 pm C0 2 channels are included in the 
temperature profile retrieval step in Version 5 because radiances in these channels are unaffected by clouds in the field of 
regard. An analogous set of temperature sounding channels is also used in Version 5.24. 


4. CHANNELS AND FUNCTION USED IN DIFFERENT STEPS OF VERSION 5.24 

Figure 2 shows a typical AIRS brightness temperature spectrum and includes the channels used in each Version 5.24 and 
Version 5.24 AO for cloud clearing, and in each of the different steps of the AIRS physical retrieval algorithm. The 
channels used in each step are summarized below. 

4.1. Temperature profile retrievals 

The location and number of channels used in the temperature profile retrieval step (red stars) are significantly different 
from those of Version 4 3 . The major difference is in the incorporation of non-LTE C0 2 temperature sounding channels 
in the spectral region 2358 cm' 1 - 2386 cm' 1 in the temperature profile retrieval step, which now uses 49 channels 
between 2197 cm' 1 and 2395 cm' 1 that are sensitive to both stratospheric and tropospheric temperatures, as well as 53 
stratospheric sounding channels between 667 cm' 1 and 713 cm' 1 . Tropospheric sounding longwave C0 2 channels are now 
used only for cloud clearing (yellow stars), and are no longer used in the temperature profile retrieval step. These same 
channels are also the ones used in the cloud parameter retrieval step 3 . Version 5.24 includes AMSU-A channels 3, 6 and 
8-14 in the temperature profile retrieval step while Version 5 AO does not. Versions 4 and 5 included AMSU-A channels 
4 and 5 in the temperature profile retrieval step, but they became noisy in 2007 and are no longer used. AMSU-A 
channel 7 was noisy at launch and was never used in any step of the retrieval process. 

4.2 Surface parameter retrievals as done in Version 4 and Version 5 

In Version 5, as in Version 4, 25 channels are used in the surface parameter retrieval step, 15 of which are between 759 
cm' 1 and 1228 cm' 1 , and the remainder are between 2450 cm' 1 and 2659 cm' 1 . These channels are used to determine 
simultaneously surface skin temperature T s , surface spectral emissivity £ v , and short-wave surface bi-directional 
reflectance p v . 
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Given an initial surface spectral emissivity guess £°, the final surface spectral emissivity e v is expressed in Version 5 as 


£ v = £ v + Z&4IKV) 


(4a) 


and the final surface bi-directional reflectance p(v) is expressed as 

Pv = Pv + Z& B t Gj(v) = ^2 + Sjl® B t C((v) . (4b) 

N F and yV G in Equations 4a and 4b are the numbers of spectral emissivity and spectral surface bi-directional reflectance 
functions being solved for in the physical retrieval step. Surface bi-directional reflectance is solved for only during the 
day. Therefore, including the surface skin temperature which is also solved for in the surface parameter retrieval step, a 
total of 1 + N F 4- N c unknown coefficients are solved for during the day and 1 + N F unknown coefficients are solved 
for at night. The emissivity perturbation functions F and G are triangles linear in frequency. In the case of a single 
function in a spectral region, the correction to the initial guess is constant in frequency. Otherwise, the correction is 
piecewise linear. 


In Version 4, N F = 2 and there is 1 spectral emissivity perturbation function solved for covering the longwave portion of 
the spectrum and 1 function covering the shortwave portion of the spectrum, with unknown coefficients A x (At w ) and 
A 2 (A sw ). In Version 5, N F = 4, with coefficients of 3 longwave emissivity perturbation functions being solved for as 
well as of 1 shortwave emissivity perturbation function. In both Version 4 and Version 5, N c = 1. 

Over land, £° is provided by a cloud cleared regression step 5 . Over non-frozen ocean, is set equal to the values found 
in the AIRS Science Team ocean emissivity model. That model is re-interpolated from calculations by Paul VanDelst 9 
using algorithms by Wu referenced in Wu and Smith 10 . 

4.3 Improved surface parameter retrieval methodology used in Version 5.24 

The basic surface parameter retrieval methodology used in Version 5.24 is very similar to that used in Version 5 with 
two exceptions. The form of the spectral emissivity and surface bi-directional reflectance perturbations has been 
modified from that shown in equation 4 to a more stable multiplicative form 


(l-e y ) = (l-e»)(l + Sil^Cr)) 

(5a) 

Pv=Pv°(l+ESi G,(v)). 

(5b) 


More significantly, 61 channels between 2395 cm' 1 and 2664 cm' 1 are now used in the surface skin temperature retrieval 
step, and 3 shortwave emissivity perturbation functions ( N F = 3) and 3 shortwave spectral bi-directional reflectance 
perturbation functions ( N G = 3) are now solved for simultaneously with T s . These channels are shown in dark blue in 
Figure 2. 79 channels, between 758 cm' 1 and 1250 cm' 1 , are now used in the subsequent surface spectral emissivity 
retrieval step, in which 3 longwave emissivity perturbation functions are solved for ( N F = 3 as in Version 5). In this 
step, T s is held fixed at the value determined from the previous skin temperature retrieval step. These channels are 
indicated by purple stars in Figure 2. 

4.4 Constituent profile retrievals 

As with earlier versions, constituent profile retrievals are performed in separate steps, each having their own set of 
channels and functions 2,3 . Figure 2 shows in different colors the Version 5.24 channels used in each of these retrieval 
steps. The q(p) retrieval (pink stars) uses 33 channels in the spectral ranges 1377 cm' 1 to 1605 cm' 1 and 2608 cm' 1 to 
2656 cm' 1 ; the 0 3 (p) retrieval (green stars) uses 41 channels between 997 cm' 1 and 1069 cm' 1 ; the CO(p) retrieval 
(green stars) uses 36 channels between 2181 cm' 1 and 2221 cm' 1 ; and the CH 4 (p) retrieval (brown stars) uses 71 
channels between 1230 cm' 1 and 1356 cm' 1 . 


6. VERSION 4 AND VERSION 5 QUALITY CONTROL 

The AIRS Science Team Version 4 retrieval algorithm introduced for the first time the concept of different quality 
control for different geophysical parameters derived from a given sounding 4 . With regard to temperature profile, there 
were 3 different flags: 1) Stratosphere Good, i.e, T(p) is good down to 200 mb; 2) Mid-Troposphere Good, i.e., T(p) is 
good down to p surf /2; 3) and Lower Troposphere Good, i.e., T(p) is good down to p surt . Version 5 further refined the 
concept so as to determine a pressure, p bes t> down to which each retrieval is considered acceptable for data assimilation 
purposes. The pressure p best is determined based on thresholds of the case-by-case error estimates of T(p), ST(p). Version 
5 uses an analogous methodology to perform Quality Control on retrieved sea surface temperatures T s , based on 
thresholds of 8 T s . Susskind et al. 5 provides more details about how error estimates are derived and used for Quality 
Control. 


7. VALIDATION OF QUALITY CONTROLLED RETRIEVALS 


7.1 Temperature profiles 

Figure 3 shows statistics for nighttime 1 km layer mean temperatures for global Quality Controlled Version 4 and 
Version 5.24 retrievals for January 25, 2003. Figure 3a shows the percent of all nighttime cases accepted in a given layer. 
Percent yields using the Version 5 Quality Control thresholds on Version 5 retrievals are shown in Figure 3a in red (dark) 
and those using the Version 4 thresholds on Version 4 retrievals are shown in gray (light). The percent accepted at 70 mb 
in Version 4 represents the percent of all cases that passed the Stratosphere Good Test. Version 5.24 accepts all cases 
above 30 mb, and applies Quality Control based on thresholds of error estimates starting at 30 mb. Results are shown 
including all cases down to p best . Quality Controlled Version 5.24 retrievals have a much higher yield beneath 200 mb 
than Version 4. Increasing spatial coverage of retrievals is very important for both data assimilation purposes and climate 
and process studies, provided their accuracy is sufficient for the appropriate application. Quality Controlled Version 5.24 
retrievals have a slightly lower yield between 100 mb and 200 mb than Version 4. This is because not all cases passing 
the Version 4 Stratospheric Temperature Test were actually good all the way down to 200 mb. Version 5.24 also has a 
higher yield above 100 mb because many cases in which the Version 4 retrieval was found not acceptable in the lower 
stratosphere were actually of high quality in the mid-upper stratosphere. 

Figure 3b shows the RMS of 1 km layer mean temperature differences between Quality Controlled Version 5.24 
retrievals (red and pink), or Version 4 retrievals (black and gray), and the collocated ECMWF analysis. Dark colors 
indicate statistics for appropriate retrieval layers selected using the Version 5 (dark) Quality Control, and the light colors 
indicate statistics for those layers selected using the Version 4 (light) Quality Control. The most meaningful comparisons 
between Version 5.24 and Version 4 retrievals are made using a common ensemble of cases, i.e., pink vs. gray or red vs. 
black. Version 5.24 retrievals are considerably more accurate than the corresponding set of Version 4 retrievals in the 
harder (dark) cases at all levels beneath 70 mb with a minor exception near 600 mb. This improvement is smaller in the 
lower yield (presumably easier cloud condition) cases accepted by Version 4 standards. 

Figures 4a and 4b show analogous results for daytime cases. There is no appreciable difference in the results shown 
during the day compared to the night. This shows that the effects of non-LTE on the radiances is well accounted for by 
the new RTA. In addition, it shows that effects of solar radiation reflected by clouds is accounted for by the cloud 
clearing methodology which produces R t . Finally, the effects of solar radiation reflected by the surface is also well 
accounted for by the retrieval of p h and it’s use in the computation of expected radiances R^clr- 

The improvement in Version 5 retrieval accuracy compared to Version 4 is a result of at least two factors, both related to 
the use of 4.2 pm C0 2 tropospheric sounding channels in the Version 5 temperature profile retrieval step compared to the 
combined use of 4.2 pm and 15 pm C0 2 channels in the Version 4 temperature profile retrieval step. The main reason for 
this choice of channels in Version 5 is that cloud clearing errors will result in smaller errors in cloud cleared brightness 
temperatures A0* for 4.2 pm channels than for the 15 pm channels. This principle allows for the determination of 
potentially more accurate temperature soundings under more difficult cloud conditions. This is demonstrated in Figures 
3b and 4b because the relative improvement in Version 5 retrieval accuracy over Version 4 retrieval accuracy is 
considerably larger in the more difficult cloud cases (dark lines) than in the easier cloud cases (light lines). This was not 
the only reason for improvement of Version 5 temperature profiles compared to Version 4 however. The inherent 
temperature sounding capability of the 4.2 pm channels may actually be better than those of the 15 pm channels, as 
indicated by the fact that Version 5 results are better in the mid-upper troposphere for easier cases as well. 

Version 5 Quality Controlled temperature profile sounding accuracy is also considerably better than that of Version 4 
between 70 mb and 200 mb. Part of this is due to the fact, alluded to previously, that not all “Stratosphere Good cases” 
are good down to 200 mb, as assumed in Version 4. These cases of poorer upper tropospheric soundings most likely 
contain very high clouds which are not appropriately screened in the Version 4 Quality Control methodology. Version 5 
retrievals between 70 mb and 200 mb using Version 4 Quality Control also perform better than those produced in 
Version 4. This appears to be partly due to better Version 5 performance under more difficult cloud conditions. In 
addition, it may also indicate better upper tropospheric sounding capability contained in the 4.2 jam sounding channels as 
compared to the 15 pm sounding channels. 


PRESSURE (mb) PRESSURE (mb) 


Percent of All Cases Included 

January 25, 2003 

G I oba I 


AM Orbits 



Version 5.24 Using V5 quality oontrol 


Figure 3a 

Percent of All Cases Included 

January 25. 2003 

Globa I 


LAYER MEAN RMS TEMPERATURE (°C) 
Differences from ECMWF 
January 25 , 2003 

G I oba I 


AM Orbits 



Version 4 Using VS quality oontrol 

Version 5.24 Using VS quality oontrol 
Version 5.24 Using V4 quality oontrol 

Figure 3b 

LAYER MEAN RMS TEMPERATURE ( 8 C) 
D i f f erences from ECMWF 
January 25 , 2003 
Global 


PM Orbits 



■ Version 5.24 Using V5 quality oontrol 


Figure 4a 


PM Orb 1 1 s 



— — Version 4 Using VS quality oontrol 

Version 5.24 Using V5 quality oontrol 

Version 5.24 Using V4 quality oontrol 

Figure 4b 


7.2 Quality Controlled ocean surface skin temperature 


The comparisons of the accuracy of Quality Controlled atmospheric temperature profiles shown in the previous section 
was made using results of Version 4 on the one hand and Version 5.24 on the other. A major difference between these 
two sets of results is due to the fact that most of the tropospheric temperature profile information in Version 4 came from 
the use of for 15 |im C0 2 channels in the retrieval process, while Version 5.24 did not use for any tropospheric 
sounding 15 Jim C0 2 channels. There were also other differences between Version 4 and Version 5.24, as well, such as 
the Quality Control methodology. The improvement in the ability to obtain accurate SST’s under more difficult cloud 
conditions by using R t for only shortwave window regions to solve for T s as compared to combined use of R t in both the 
longwave and shortwave window regions, is best isolated by a comparison of results of Version 5.24 and 5.0, which are 
otherwise identical to first order. 

Figure 5 shows histograms of the counts of errors of quality controlled ocean SST retrievals as a function of SST 
differences from ECMWF “truth” on January 25, 2003, obtained using Version 5, and Version 5.24. The mean difference 
from ECMWF, the spatial standard deviation of the difference, the percent of all cases accepted, and the percent outliers 
(errors more than 3K from the mean) are indicated in Figure 5. Results for day and night are shown separately. Version 
5.24 Quality Controlled SST retrievals have considerably higher yields with similar standard deviation and percent 
outliers than the comparable Quality Controlled Version 5 SST’s. The improvement is actually greater during the day 
than at night. 
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Figure 6 shows the percent of all SST cases 50°N-50°S accepted as a function of retrieved fractional cloud cover for both 
Version 5 and Version 5.24. Results are shown separately for day and night cases. Version 5 Quality Control 
methodology automatically rejected all retrievals with fractional cloud cover greater than 90%, both with regard to 
temperature profile and surface skin temperature. Version 5.24 does not have this restriction. All cases shown in this 
figure were included in the statistics shown in Figure 5. As shown in Figure 5, Version 5.24 has a considerably higher 


percent yield than Version 5, with comparable overall accuracy, both under daytime and nighttime conditions. Version 
5.24 Quality Controlled SST retrievals have higher yields than Version 5 at all cloud fractions both day and night. The 
relative improvement in the ability to determine accurate Quality Controlled SST’s between Version 5.24 and Version 5 
increases dramatically with increasing cloud cover. For example, at 50% fractional cloud cover, Version 5.24 accepts 
roughly 34% of the cases during the day and 16% at night, while Version 5 accepts roughly 8% of the cases under both 
day and night conditions. This increased yield under more difficult cloud conditions is significant for climate 
applications. 
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8. SUMMARY 

This paper shows that improved AIRS Quality Controlled soundings of atmospheric temperature profile T(p) and ocean 
surface skin temperature T s are obtained by use of cloud cleared radiances of only shortwave window channels in the 
determination of T s , and primarily of 4.2 jim CO 2 channels in the determination of T(p). Observed radiances R t are used 
in the longwave 15 jim C0 2 band and lljim window in the determination of for all channels. This improvement in 
sounding capability exists for relatively clear situations, and is magnified with increasing cloud cover. Results obtained 
during the day are at least as good, if not better, than those obtained at night. Similar results are obtained in the Version 
5.24 AIRS Only system (not shown in this paper). 

This improved sounding methodology is made possible as a result of a few factors, both theoretical and instrumental. 
First of all, use of improved radiative transfer physics allows for the accurate computation of the effects of non-LTE on 
the observed radiances in the 4.2 jim C0 2 band during the day. Secondly, the cloud clearing methodology used to 
analyze AIRS data removes the effects of solar radiation reflected by clouds from the clear column radiances used to 
generate the solution. Effects of solar radiation reflected by the surface on the shortwave radiances are accounted for 
directly as part of the surface parameter retrieval process. Also, improved error estimates and Quality Control 
methodology allows for the screening of cases in which aspects of the retrieval methodology may have performed 
poorly. 


Finally, and most significantly, use of these theoretical improvements is made practical because the noise in the AIRS 
shortwave channels is extremely low. This approach is not practical using IASI data because the IASI noise is too high in 
the shortwave portion of the spectrum. The findings of this paper are significant with regard to design considerations of 
future high spectral resolution IR sounders, especially for geostationary (GEO) orbit. It is optimal on such sounders to 
have spectral coverage beyond 2400 cm' 1 with very low noise, to produce accurate soundings under most cloud 
conditions. It is not essential for advanced IR GEO sounder to be accompanied by a GEO microwave sounder to achieve 
these results. 
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